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Q"^ , find that the black hole mass scales at threshold like Mbh oc \p — p*\'^ , "where 
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*Electronic address: yoav@csd.uwm.edu 
^Electronic address: hose@csd.uwm.edu 
■'"Electronic address: leonard@cosmos.phys.uwm.edu 

1 



The discovery of universality and scaling at the onset of black hole formation [|l| may have 
important implications in understanding the structure of solution-space in gravitational the- 
ory. Choptuik studied numerically the collapse of a spherically symmetric self-gravitating 
real scalar field in four- dimensional (4D) Einstein gravity and considered one-parameter 
families of initial data, S[p], where p is a parameter specifying the strength of the gravita- 
tional self- interaction of the scalar field. He found that there exists a critical value, p = p*, 
such that for p < p*, the "subcritical case", no black hole is formed and the solution is 
regular, while for p > p*, the "supercritical black hole is formed. Furthermore, by 

careful analysis of the solutions near criticality, p = p*, he found that as p approaches p* 
from above, the mass of the created black hole approaches zero, and when a black hole just 
forms, its mass scales as Mbh oc \p — p*\'^, where the critical exponent 7 ^ 0.37. The critical 
solutions also exhibit discrete self-similarity Similar behavior was found in other models 
of non-linear gravity 0]. 

In the semiclassical scenario, i.e., for a quantum field propagating on a classical dynamical 
background metric, the created black hole of mass M radiates in 4D with the Hawking 
temperature Th oc M'^ §. As M ^ near criticality, Th becomes large and quantum 
effects are clearly important. To gain some understanding of how quantum effects may alter 
the classical scenario near criticality, we investigate a tractable two-dimensional (2D) model 
that exhibits classical scaling as well as significant quantum effects near criticality. 

The model we study is 2D dilaton gravity coupled to a massless scalar field. We first 
consider the classical theory and then include quantum effects. We find that classically 
there is universal scaling of the black hole mass near criticality, Mbh oc \p — p*\'^, where 
7 ~ 0.53 ± 0.01 is independent of initial data. In addition, we find that the ground state of 
the classical theory gives a lower bound on the energy of spacetimes resulting from dynamical 
processes in which a black hole is formed. This is related to a radiation energy deficit, Arad, 
that we describe later. At the semiclassical level, although the Hawking temperature in our 
2D model is independent of M, we find that quantum effects alter significantly the critical 
behavior. Most interestingly, when p is sufficiently close to p*, the classical scaling behavior 
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breaks down and Mbh approaches a non-zero constant value that depends on the initial 
data. The existence of this mass gap at the quantum level is a new result that may be of 
fundamental interest. 

The classical theory is described by the 2D Callan-Giddings-Harvey-Strominger (CGHS) 
action Q] 

-^E(v/.)^|, (1) 

1=1 

where is the dilaton field, TZ is the 2D Ricci scalar formed from the metric tensor Q/j^i^, A is a 
positive constant, and the fi are massless scalar matter fields conformally coupled to the 

2D geometry. We work in the conformal gauge, g = 5'++ = and g^^ = — (1/2) exp(2p), 

(i.e., ds"^ = — exp[2p{x~^ , x~)]dx~^ dx~) , where {x~^,x~) are the "Kruskal" null coordinates in 
which (j){x^,x^) = p{x^,x^) p. 

The general vacuum solution of the classical theory (|l]) is = —(1/2) ln(— A^a;"'"x~ + C), 
where AC is its ADM mass [0. For C > the vacuum solutions describe static 2D black 
holes. The C = solution is the linear dilaton vacuum (LDV), which is the ground state 
of the theory. Solutions with C < have timelike naked singularities in the strong coupling 
region where exp(20) — > 00. The weak coupling, asymptotically fiat region will be taken to 
be on the right-hand-side (rhs) of the spacetime [^. In this work we take C < 0, but we avoid 
the region of strong coupling and the singularity by considering the solutions only in the 
weak coupling region exp(20) < exp(20c) for a given constant (pc, and by imposing reflecting 
boundary conditions on the timelike boundary curve (f){x^,x^) = (pc- A previous study of 
such a system at the onset of black hole formation was based on boundary conditions 
that mix classical and one-loop contributions and do not have a standard classical limit. 
The reflecting boundary conditions that we impose here have a standard classical limit. 
Let the boundary be described by the curve x~^ = x~^{x~). Then the reflecting boundary 
condition is 

2^-(^-)= (2) 
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'S'cGHS = TT- I ^^2; A 



,-24, 



7^ + 4(V0)2 + 4A2 



where T±± = (1/2) J2iid±fiY are the components, in the {x~^,x~) coordinates, of the stress 
tensor of the massless scalar fields. The general solution for the conformally coupled matter 
fields is fi{x~^, x~) = + The initial data at right asymptotic past null infinity, 

QJi, are therefore given by fi~{x~^) or T++(a:;+). Working in the large limit [^] we are able 
to choose the boundary curve to be at exp(— 20^) 0^ and get a second order non-linear 
ordinary differential equation (ODE) for the classical boundary curve 

x-+P+{xt)/X' 

where ^4.(2;+) = J^~^ T^^{x^)dx^ . We solve the ODE (§) numerically for different initial 
data, T++(a;^), with compact support, xf < x"^ < a;^. To the past of xf, i.e., for x+ < xf, 
we have a vacuum solution = —(1/2) ln(— A^x^x" +C), with C being a negative constant. 
In this region we have an analytical solution for the boundary curve: x^(x^) = C/x^ , where 
x" < C/xi . For x~ > C/xi, we integrate (|]) numerically to find the boundary curve and 
the corresponding solutions for the metric and dilaton field. 

We study in detail the two families of initial data shown in Fig. |1|. Family (a) is a 
two-parameter family specified by the amplitude. A, and the inverse width, F, of the profile 
T^y{v), where v = X^^ In(Ax^) and u = ln(— Ax^ — P_|_(oo)/A) are the manifestly asymp- 
totically flat null coordinates on 53^. Family (/3) is a five-parameter family described by Ai, 
A2, Fi, F2 and 6. All the parameters are scaled by appropriate powers of A to make them 
dimensionless. In these two cases (and in other cases that we studied) we find that fixing 
all but one arbitrarily chosen parameter, say p, yields a regular evolution with no black hole 
formation if and only if p < p*. The critical value, p*, depends on the values of the other 
fixed parameters. As in the 4D solution-space, the only "intermediate solution" separating 
black hole solutions from regular ones is the critical solution with p = p*. 

Next we show that there is a universal scaling of the black hole mass near criticality in 
the classical theory. In order to find the mass of the created black hole we calculate the 
outgoing radiation reaching right asymptotic future null infinity, ^3^, after being reflected 
from the boundary. This radiation is described by the stress tensor, T„„(-u). Then we find 
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that the Bondi mass at late times {u —>■ oo) on is M^^^^^ = MBondi(^^ —>■ oo) = Madm — 



Tuu{u)du. Here Madm is the ADM mass ||T0[ of the spacetime, defined such that the LDV 



ground state has zero ADM mass. In Fig. |^ we show the numerical results for the reflected 
outgoing radiation, Tuu{u), for two nearby sets of initial data. The two solutions correspond 
to family (a) with the same width, F = 0.2, but with different amplitudes, A. The upper gray 
curve corresponds to a subcritical solution near criticality, A = As^h = A* — ei = 1.546. The 
black curve corresponds to a supercritical case near criticality, A = Agupcr = A* + €2 = 1.547, 
where A* is the critical value of the amplitude parameter A. 

Although the initial data are very similar in the two cases, the late-time outgoing ra- 
diation is quite different. For the subcritical solution near criticality, the boundary curve 
becomes almost null and the late-time outgoing radiation is highly blue-shifted In the 
supercritical case on the other hand, this blue-shifted radiation is absent, and we flnd that 
part of the incoming radiation does not reach future asymptotic null inflnity. Since Tuu is 
quite different in the two cases, the Bondi mass deflned through the integral jTuu{u)du is 
obviously different. For all the subcritical solutions the Bondi mass at late times, Mg'^jj, is 
equal to the negative ADM mass of the initial spacetime, AC. However, for the supercritical 
solutions Mg'^j; is always positive. We deflne the "radiation energy deflcit", Arad, to be 
Arad = hme_^o {m^I^^[p = p* + e] - m'^I^^Ip = p* - e]|. The fact that A^ad turns out to be 
non-zero is a striking manifestation of the non-linearity of our problem. 

We find that Arad does not depend on the profile of the infalling matter, but only on the 
initial vacuum geometry. Explicitly, 

Arad = -AC = A|C| , (4) 

where C is the constant specifying the initial vacuum geometry. The existence of this non- 
zero Arad can be interpreted as implying that in the classical theory the LDV is a "stable 
ground state" . In particular, we have started with a negative-mass geometry and find that 
the process of throwing in matter to form a black hole results in spacetimes having non- 
negative mass relative to the LDV. In agreement with cosmic censorship, this suggests that 



systems in this classical theory will not evolve to states of energy lower than the LDV, which 
have naked singularities. Once the infalling matter configuration is dense enough, p > p*, 
part of the incoming energy is "swallowed" by the strong curvature region to compensate 
for the negative mass of the initial spacetime and make the future spacetime a solution 
with positive energy compared to the LDV. For our initial geometry with negative ADM 
mass AC, the minimum energy that must be swallowed in order to get a positive-energy 
spacetime, is precisely Arad- The resulting positive-energy spacetime is a 2D black hole. We 
find that in this classical theory the mass of the created black hole, Mbh = M^^^^^lp > p*], 
approaches zero as p approaches p* from above. Moreover, the black hole mass scales as 

MbhOc|p-p*r (5) 

at threshold, where 7 = 0.53 ± 0.01. The value of 7 does not appear to depend on the 
profile of the infalling matter or the initial geometry. In Fig 3(a) we plot ln(Mbh/A) + 
versus ln\pi — p*\ for different parameters pi. The constants are chosen such that the 
three curves intersect at a given point. The slopes of all the lines are the same within our 
numerical accuracy. 

Next we consider quantum effects. We quantize the scalar fields on a classical dynam- 
ical background geometry consisting of the metric and the dilaton field. The semiclassical 



effective action that we study is |11 



'S'cfT = 5'cGHS ~ J d^x^-g{x) j (fx'^-g{x')n{x)G{x, x')lZ{x') 



+ ^ / d\ [(V0)2 - 07^] , (6) 

where k = Nh/12 and G{x,x') is an appropriate Green function for the massless scalar 
fields. The second term on the rhs of Eq. is the Polyakov-Liouville effective action 
derived from the trace anomaly of the 2D massless scalar fields |]12[, and the last term on 



the rhs of Eq. is a local counter-term that we add to our effective theory in order 
to make it exactly solvable |TT[|. The vacuum solutions of our semiclassical theory are 



— (1/2) ln[— A^x'^x — (fi;/4) ln(— A^x"''x ) + C], where C is a constant. The ground 
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state is the solution with C = Cq = (K/4)[ln(K/4) — 1]. As in the classical case, we impose 
reflecting boundary conditions on a timelike boundary curve and study the solutions only 
in the weak coupling region. The reflecting boundary condition is a modiflcation of Eq. 
where T±± are now components of the total stress tensor, including the classical and one- 
loop contributions. There is also an extra term due to possible particle creation from the 
boundary, which is effectively a moving mirror One can eliminate the moving mirror 

term by taking the large N limit together with exp(— 20c) 0+ g. In contrast to (|), the 
resulting nonlinear ODE for the semiclassical boundary curve, = is 

2 



A2 AX^x^ 



_|_ o *^*^-^ _|_ 

dx~'^ dx~ 



doc ^ 
dx~ 



0. (7) 



We consider initial data with no quantum radiation on and classical incoming radi- 
ation corresponding to the proflles shown in Fig. |l|. We flnd that as in the classical case, 
the solutions are regular and no black hole is formed if and only if p < p* (where p* is the 
critical value of the parameter p in the semiclassical case). For p > p*, the black hole that 



is formed evaporates by emitting Hawking radiation W^. To look for scaling analogous to 
Eq. (^, we calculate the initial mass of the created black hole. It is this mass that reduces 
in the limit Nh to the mass of the classical black hole that appears in Eq. (^. The 
initial black hole mass, mS, is the Bondi mass at M = where uq is the minimum of the 
apparent horizon u = nah(^)- The apparent horizon is the solution of the equation d+(f) = 0. 
Explicitly, M^^ = Madm — /"J^ T^udu, where T„„('u) is the total stress tensor of the outgoing 
radiation and Madm is the ADM mass of the spacetime, deflned such that the semiclassical 
ground state, i.e., the static vacuum solution with C = Cq, has zero ADM mass. 

We examine quantum effects by considering cases with different values of k = Nh/12. 
In Fig. 1(b) we plot \xi{M^l/X) versus In |p — p*| for two different values of /t. We do so in 
the case of the family {a) of initial data shown in Fig. [1|, where the free parameter p is the 
amplitude A, and the parameter F = 0.2 is flxed. For large values of M^^ the initial mass of 
the black hole behaves like that of the classical black hole (|^). However as M^^^ decreases, 
deviations from the classical behavior appear. Unlike the classical case, as p approaches 



p* the initial mass, M^^ of the created black hole does not generally approach zero, but 
rather approaches a constant, Mgap. We find that this mass gap, Mgap, depends not only 
on the value of n but also on the initial data. For the cases shown in Fig. ^b) we have 
Mgap/A ~ 0.014 for k = lO'^ and Mgap/A ~ 0.0037 for k = IQ-^. The radiation energy 
deficit in the quantum case is A^^^^ = \{Co ^ C) + Mgap. 

The deviations from the classical scaling law (|^) become significant for values of M^^ that 
are of order kX. When M^^^ takes values that are of order or less than nX, the semiclassical 
approximation remains valid provided is sufficiently large However for any fixed 



value of A^, no matter how large, the semiclassical approximation breaks down when p is 
sufficiently close to p* P|, i.e., as one moves to the far left well beyond the region shown in 
Fig. §(b). Thus, determining the behavior of the curve as ln|p — p*| — oo requires full 
quantization of the theory. We are currently studying the phase structure of the semiclassical 
theory near threshold which appears to be richer than that in the classical case. 

In this work we investigate numerically the collapse of massless scalar fields in 2D dilaton 
gravity. We find that classically the black hole mass at threshold obeys a power-law, Mbh oc 
\p — p*\'^, where 7 = 0.53 ± 0.01. However, quantum effects give rise to a mass gap that 
depends on the initial data. 

We would like to thank B. Allen, J. Friedman, D. Garfinkle, T. Piran, and A. Steif 
for very helpful discussions. This work was supported by the National Science Foundation 
under Grant No. PHY 95-07740. 
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FIGURES 



FIG. 1. Two families of initial data described by the stress tensor T^„. 

FIG. 2. The stress tensor T^u, in units of A^, describing the outgoing radiation. The upper gray 
curve corresponds to a subcritical solution just below criticality and the black curve corresponds 
to a supercritical solution just above criticality. 

FIG. 3. ln(Mbh/A) versus \n\p — p*\ in the classical case, (a), and the quantum case, (b). 
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